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Precessional instability in binary black holes with aligned spins 
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Binary black holes on quasicircular orbits with spins aligned with their orbital angular momentum 
have been test beds for analytic and numerical relativity for decades, not least because symmetry 
ensures that such configurations are equilibrium solutions to the spin-precession equations. In this 
work, we show that these solutions can be unstable when the spin of the higher-mass black hole is 
aligned with the orbital angular momentum and the spin of the lower-mass black hole is antialigned. 
Spins in these configurations are unstable to precession to large misalignment when the binary 
separation r is between the values r u a± = (y^xT ± ^/<7X'2) 4 (1 — g) _2 M, where M is the total mass, 
q = vrii/mi is the mass ratio, and xi (X 2 ) is the dimensionless spin of the more (less) massive black 
hole. This instability exists for a wide range of spin magnitudes and mass ratios and can occur in 
the strong-field regime near the merger. We describe the origin and nature of the instability using 
recently developed analytical techniques to characterize fully generic spin precession. This instability 
provides a channel to circumvent astrophysical spin alignment at large binary separations, allowing 
significant spin precession prior to merger affecting both gravitational-wave and electromagnetic 
signatures of stellar-mass and supermassive binary black holes. 

PACS numbers: 04.25.dg, 04.70.Bw, 04.30.-w 


Introduction. - Black holes (BHs) have been ob¬ 
served in two distinct regimes: stellar-mass BHs 
(5M 0 < to < 100-Mq) accrete from companions in x-ray 
binaries m, while supermassive BHs shine as quasars 
or active galactic nuclei (AGN) j3J j3l. Both types of 
BHs naturally occur in binaries: the massive stellar pro¬ 
genitors of stellar-mass BHs are typically formed in bi¬ 
naries, while supermassive BHs form binaries following 
the mergers of their host galaxies [B]. Gravitational ra¬ 
diation circularizes the orbits of these binaries [7; and 
causes them to inspiral and eventually merge, making 
them promising sources of gravitational waves (GWs) for 
current and future GW detectors I8HT5|. The spins of 
these binary BHs need not be aligned with their orbital 
angular momentum: stellar-mass BHs may recoil during 
asymmetric collapses tilting their spins with respect to 
the orbital plane [TBl - fTH] . while the initial orbital plane 
of supermassive BH binaries reflects that of their host 
galaxies and is thus independent of their spin. Gravi¬ 
tational effects alone will not align the BH spins with 
the orbital angular momentum mm, but astrophys¬ 
ical mechanisms exist that drive the BH spins towards 
alignment in both regimes. The first BH to collapse in 
stellar-mass BH binaries may accrete in a disk from its 
as yet uncollapsed companion, while both members of 
a supermassive BH binary may accrete from a common 


circumbinary disk. Warps in these accretion disks can 
align the BH spins with the orbital angular momentum 
[2fH23| . but if the initial misalignment between the BH 
spin and accretion disk is greater than 90°, the BH may 
instead be driven into antialignment [23] . 

Misaligned spins cause the orbital angular momentum 
to precess j25H2T[ . modulating the emitted GWs (25] , 
Spin misalignment is both a blessing and a curse for GW 
data analysis: it increases the parameter space of tem¬ 
plates needed to detect GWs via matched filtering but 
also breaks degeneracies between estimated parameters 
in detected events (25j. Misaligned spins at merger can 
generate large gravitational recoils [30H32] . ejecting su¬ 
permassive BHs from their host galaxies. Spin precession 
may also be responsible for the observed X-shaped mor¬ 
phology of AGN radio lobes [331 S3 • Given the impor¬ 
tance of spin misalignment, it is worth investigating the 
robustness of aligned spin configurations. In the general 
case that the BHs have unequal masses, there are four 
distinct (anti-)aligned configurations, which we refer to 
as up-up, up-down, down-up, and down-down. The direc¬ 
tion before (after) the hyphen describes the more (less) 
massive BH and up (down) implies (anti-) alignment of 
the spin with the orbital angular momentum. By sym¬ 
metry, all four configurations are equilibrium solutions to 
the orbit-averaged spin-precession equations m, but are 
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FIG. 1. Effective-potential loops £±(S) for binary BHs with mass ratio q = 0.9, dimensionless spins xi = 1, X 2 = 0.14, and total 
angular momentum J = \L + Si — S 2 I, corresponding to the up-down configuration. For binary separations r > r u[ j+ — 337 M 
(left panel), the up-down configuration at 5 m i n marked by a red circle is also a minimum (marked by the lower triangle). At 
intermediate separations r u d+ > r > r u< j_ — 17 M (middle panels), misaligned binaries with the same value of the conserved £ 
exist along the dashed red line. Perturbations SJ, S £ will cause S to oscillate between the points S± where this line intersects 
the loop, making the up-down configuration unstable. For r < r u d- (right panel), the up-down configuration is again a stable 
extremum, now a maximum (marked by the upper triangle). An animated version of this figure is available online at Ref. [35) . 


these solutions stable? To answer this question, we inves¬ 
tigate how the configurations respond to perturbations of 
the spin directions using our recently developed approach 
for studying generically precessing systems (MUSH- As 
we will demonstrate below, the up-down configuration is 
unstable for certain choices of binary parameters, with 
significant consequences for GW data analysis and astro¬ 
physics. 

Generic spin precession. - Here we briefly summa¬ 
rize the approach to spin precession described in de¬ 
tail in [Ml 133 using units where G = c = 1. Bi¬ 
nary BHs with total mass M = toi + m 2 , mass ratio 
q = m 2 /mi < 1, symmetric mass ratio rj = q/(l + q) 2 , 
and spins S, = evolve on three distinct time 

scales: the orbital time t or b = (r 3 /M) 1 / 2 on which 
their separation r changes direction, the precession time 
tpre = (t m -b/ri){r/M) on which the spins and orbital an¬ 
gular momentum L change direction, and the radiation- 
reaction time £rr = {t OI b/v)( r /AT) 5 / 2 on which the mag¬ 
nitudes r and L decrease due to GW emission. The 
relative orientations of the spins are often specified by 
the two angles cos &i = Si • L and the angle A<f> be¬ 
tween the projections of the two spins onto the orbital 
plane, all of which vary on t pre . The spin orientations 
can equivalently be specified by the magnitudes of the 
total spin S = Si + S 2 , the total angular momentum 
J = L + S, and the projected effective spin 2H 132] 
£ = M _2 [(l+g)Si-t-(l- 1 -g -1 )S 2 ]-L. This specification has 
the advantage that only S evolves on f pre , while J evolves 
on Irr and £ is conserved throughout the post-Newtonian 
(PN) stage of the inspiral (r > 10M) by orbit-averaged 
2PN spin precession and 2.5PN radiation reaction uni¬ 
on the precession time, the spin magnitude S simply os¬ 


cillates back and forth between the two roots S± of the 
equation £ = £±(S), where 

£ ± (S) = {{J 2 -L 2 - S 2 )[S 2 ( 1 + q) 2 - (S 2 - Sf)(l - q 2 )} 

± (1 - q 2 W[J 2 -(L-S) 2 }[(L + S) 2 - J 2 } 

x V[S 2 - (Si - S 2 ) 2 ][(S 1 + S 2 ) 2 - S 2 ]}/(AqM 2 S 2 L ), 

(1) 

are the effective potentials for BH spin precession. Note 
that S is the only quantity on the right-hand side of 
Eq. Q changing on f pre ; in the absence of radiation re¬ 
action, the spins return to their initial relative orientation 
after a time r(L, J, £) during which L, Si, and S 2 pre- 
cess about J by an angle a(L, J, £). The two potentials 
£±(S) form a closed loop in the 5£ plane, implying that 
the two roots S± coincide at the extrema £ m in,max(A, J) 
of the loop. At these extrema, also known as spin-orbit 
resonances m, S does not oscillate and L, S 1; and S 2 
all remain coplanar on the precession time. 

Stability of aligned configurations. - We begin with 
the up-up and down-down configurations, for which 
J = |L ± (Si + iS 2 )|, respectively. According to Eq. 0, 
the effective-potential loop reduces to a single point in 
this limit which is necessarily an extremum: S cannot 
oscillate consistent with conservation of J and £. Now 
consider the down-up (up-down) configurations for which 
J = \L — Si + S 2 I (J = |L + Si — S 2 I). The effective- 
potential loop £±(S) encloses a nonzero area for these 
values of J, implying that oscillations in S are possi¬ 
ble, except at the extrema £ m in,max- Since the spins are 
antialigned with each other in both configurations, S is 
minimized at S m i n = |Si — S 2 | and both configurations 
sit on the leftmost point of the loop, where £+(S) and 
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FIG. 2. The angles cos 6i = Si ■ L for spin-orbit resonances 
[extrema of £±(S)j for BHs with q = 0.95, Xi = 0-3, and 
X 2 = 1. The solid (dashed) curves indicate the A4> = 0 (7r) 
family and the five curves for each family correspond to binary 
separations r/M = 3000,720,170,40, and 10. The up-down 
configuration (bottom right corner) belongs to the A<f? = 0 
family for r > r u a+ — 2149M, to the A$ = 7r family for r < 
r u d- — 13M, and is unstable for intermediate values r u d- < 
r < r u d+. An animated version of this figure is available 
online at Ref. [35j . 


FIG. 3. Precession-averaged radiation reaction dJ/dL as a 
function of J and £ for binaries with q = 0.8, xi = \2 = 1, and 
separation r = 10 M in the unstable region r u d- < r < r u d+. 
Spin-orbit resonances including the up-up, down-down, and 
down-up configurations are extrema of £±{S) and constitute 
the boundary of the allowed region. All four aligned con¬ 
figurations are maxima where d,J/dL = 1, but the unstable 
up-down configuration (shown in the inset) is a cusp. An 
animated version of this figure is available online at Ref. ESI- 


£~(S) coincide. Whether this point is also an extremum 
£min,max depends on the slopes of these two functions 
at that point. Both slopes are always negative for the 
down-up configuration, implying that it is a maximum 
£max and thus a spin-orbit resonance like the up-up and 
down-down configurations. At large binary separations 
r, the slopes of £±(S) are both positive for the up-down 
configuration, making it a minimum £ min . However, be¬ 
low r ud + given by 


r u d± 


(Vxl ± Vox 2) 4 T 

(1 - q ) 2 


( 2 ) 


the slope of £_(>!?) becomes negative and up-down is no 
longer an extremum of the effective-potential loop, as 
seen in Fig. [l] At separations below r m j_, the slope of 
£,+ {S) also becomes negative and up-down is again an 
extremum, this time a maximum £ max . Misaligned BHs 
with the same values of J and £ as the up-down con¬ 
figuration but S > iSmin exist in the intermediate range 
r ud - < r < r ud +, as shown by the dashed red line. These 
misaligned BHs have an infinite precessional period r: 
they exponentially approach the up-down configuration 
on the precession time t pre but never reach it. 

The evolving relationship between the up-down con¬ 
figuration and the spin-orbit resonances parameterized 
by the angles is seen in Fig. [5] The solid curves 
show the A<F = 0 resonances [£ m i n (</)] for separations 
10M < r < 3000 M, while the dashed curves show the 
Ad> = 7r resonances [£ m ax (</)]■ The up-down configura¬ 


tion is located in the bottom right corner of this figure. 
For r > r u d+, the up-down configuration lies on the solid 
curves and belongs to the A<F = 0 family, but for smaller 
separations these curves detach from the bottom right 
corner, and thus up-down is no longer a minimum of 
£±(S'). The dashed curves indicating the Ad* = ir fam¬ 
ily migrate to the right with decreasing separation and 
reach the bottom right corner, making the up-down con¬ 
figuration a maximum of ( S ), for r < r ud - • The up-up 
and down-down configurations (top right and bottom left 
corners) belong to both resonant families, reflecting the 
degeneracy of the effective-potential loop as a single point 
that is both minimum and maximum. The down-up con¬ 
figuration (top left) always belongs to the A$ = ir family 
and is thus a maximum £ max . 

The stability of a system is determined by its re¬ 
sponse to perturbations, in this case to the spin angles 
(<5^i, <5^2, <5A<3>) or equivalently to the angular momenta 
(i 5S,SJ,S^). After such a perturbation, configurations 
that are extrema of £,±{S) (all aligned configurations ex¬ 
cept up-down for r ud - < r < r u d+) will undergo oscil¬ 
lations in S (and thus the three spin angles) that are 
linear in the perturbation amplitude, and have a period 
r that is independent of this amplitude. This is a sta¬ 
ble response equivalent to that of a simple harmonic os¬ 
cillator. The response of the up-down configuration for 
r u d- < r < r u d+ is very different, as seen in the middle 
panels of Fig. [I] S oscillates between the turning points 
S± independent of the perturbation amplitude, but the 
period r of these oscillations - as predicted by Eq. (27) of 
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Ref. diverges logarithmically as this amplitude ap¬ 

proaches zero. This is an unstable response: the time it 
takes for a zero-energy particle with dx/dt < 0 to travel 
from finite Xq to Sx in the unstable potential V = —\kx 2 
similarly diverges logarithmically with Sx. 

A perturbative analysis of nearly aligned configura¬ 
tions |41| can identify that perturbations can oscillate 
at complex frequencies (indicating an instability) in the 
same region r u < j_ < r < r ud _|_ found here, but such analy¬ 
sis cannot predict the amplitude of these perturbations or 
their response to precession-averaged radiation reaction. 
Radiation reaction. - We have shown that for 
r uc j_ < r < r ud+ , spin configurations with J and £ in¬ 
finitesimally close to the up-down configuration can expe¬ 
rience finite-amplitude oscillations in S and the angles B \, 
02, and A4>. We now investigate how these configurations 
evolve on the longer radiation-reaction time £rr. Since £ 
is conserved throughout the inspiral and L monotonically 
decreases at 2.5PN order, the only challenge is to evolve 
J. In Refs. [2(2 07j we derived a precession-averaged ex¬ 
pression for dJ/dL , a contour plot of which is shown in 
Fig-i The shaded region shows the allowed values of J 
and £ for this mass ratio, spin magnitudes, and binary 
separation. The spin-orbit resonances, being extrema of 
£±(S), constitute the boundaries of this region. The up- 
up, down-down, and down-up configurations, being spin- 
orbit resonances, lie on these boundaries. At r ud _|_, the 
up-down configuration detaches from the right bound¬ 
ary of this region [it stops being a minimum of £±(5)] 
and begins to migrate leftwards through the allowed re¬ 
gion, eventually reattaching to the left boundary at r u( j_ 
[where it becomes a maximum of (S')]. This is just an 
alternative visualization of the four panels of Fig. [lj 

For all four aligned configurations, J and L are aligned 
so dJ/dL = 1 is maximized. However, the nature of 
these maxima is very different for the stable and un¬ 
stable configurations. For the stable configurations, the 
partial derivatives of dJ/dL with respect to J and £ re¬ 
main finite, implying that neighboring points separated 
by (SJ, St;) slowly drift away at a rate that scales linearly 
with these infinitesimal quantities. The unstable config¬ 
uration however is a cusp where these partial derivatives 
approach ±oo, depending on whether this point in the 
plane is approached from below or above. Neighboring 
points (experiencing large-amplitude oscillations in S, as 
seen in the middle panels of Fig. [TJ rapidly deviate from 
the up-down configuration as it sweeps across the allowed 
region. This is an essential point: even if the stability of 
the up-down configuration is restored in the PN regime 
(r u d- > 10M), radiation reaction during the inspiral be¬ 
tween r u( j± will drive BHs initially in this configuration 
to large misalignments prior to merger. The migration 
of the up-down configuration through the plane also 
reconciles the instability with the empirical result that 
isotropic spin distributions remain isotropic during the 
inspiral [HI1D]: although nearby binaries may indeed 



FIG. 4. Normalized GW Fourier amplitude h (cf. Ref. [41] 1 
as a function of orbital frequency / and binary separation r 
during the inspiral of BHs with q = 0.75 and xi = X 2 = 0.9. 
At the initial separation r = 1000A7, the spins are nearly 
in the up-down configuration, but this configuration becomes 
unstable below r u d+ — 157M, after which large precession- 
induced modulations occur at frequencies accessible to GW 
detectors. 

be left behind, the unstable configuration will always en¬ 
counter a fresh supply, until it is restored to stability at 
the left edge of the allowed region. 

GW astronomy. - Binaries with separations in the un¬ 
stable region between r uc j± emit GWs with frequencies in 
the range / ud ± ~ 6.4 x 1O 4 Hz(M/M 0 ) _1 (1 - g) 3 /( 1 /xl± 
\fdX 2 ) 6 ) within or below the sensitivity band of existing 
and planned GW detectors In Fig.|4j we show the 

waveform of one such binary initially near the up-down 
configuration before entering the unstable region. Once 
the binary crosses the threshold at r ud +, its waveform 
develops large-amplitude precessional modulation on the 
precession time f pre . The amplitude of this modulation 
is independent of the initial deviation from the up-down 
configuration: it is set by the finite-amplitude oscilla¬ 
tions in S seen in the middle panels of Fig. |T| Modula¬ 
tion occurs on two distinct time scales associated with 
the precession of L in a frame aligned with J. In this 
frame the direction of L is specified by the polar angle 
cos 9l = L • J and the azimuthal angle 4 >r in the plane 
perpendicular to J. The longer of these time scales is r 
(the period of oscillations in 0 r), while the shorter time 
scale is (27r/a)r (the precession-averaged time for to 
change by 27r) [55(157]. Measuring this modulation could 
yield insights into the astrophysical origins of binary BHs 
m i33- Spin precession could also affect the electro¬ 
magnetic counterparts to BH mergers [42] [43] and the 
probability of ejecting a supermassive BH from its host 
galaxy [50H5121 Hil . We look forward to confronting these 
predictions with observations in the dawning age of GW 
astronomy. 
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